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Abstract 

We formulate the many-body spin dynamics at high temperature within 
the non-equilibrium Keldysh formalism. For the simplest XY interaction, 
analytical expressions in terms of the one particle solutions are obtained 
for linear and ring configurations. For small rings of even spin number, 
the group velocities of excitations depend on the parity of the total spin 
projection. This should enable a dynamical filtering of spin projections 
with a given parity i.e. a spin projection chromatography. 
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1 Introduction 

The evolution of a local spin excitation in a system of interacting spins re- 
mains a problem of broad interest in theoretical and experimental physics. On 
one extreme, it seeks the quantification of the hydrodynamic limit of spin "dif- 
fusion" PIE]; on the other, the growing field of quantum information processing 
needs a dynamical[SJ |S] control of individual spins in systems of intermediate 
size pushing the theory to its quantum limit. In spite of the naive expectation 
that quantum dynamical interferences might be limited to simple systems close 
to its ground state, most of the advances were originated on Nuclear Magnetic 
Resonance (NMR) [7] which deals with spin-spin interactions much weaker than 
the thermal energy ksT. The discreteness of the spin system produces quantum 
beats that allow an efficient cross-polarization among different spin species 
In small bounded spin systems, such as molecules, the initial spreading out of 
a local excitation is followed by a revival known as Mesoscopic Echo^\ (ME) 
manifesting the reflection at the molecular boundaries. This echo was observed 
in the proton system of ferrocene[Tfj]. (C5H 5 ) 2 Fe, and benzene^J, CqRq, ori- 
ented molecules. Indeed, the complex interferences produced by the many-body 
dipolar interactions conspire against the spin wave behavior and produce a rapid 
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attenuation of the ME. Therefore, a clean observation of the ME requires a re- 
duction of the complex Heisenberg or Dipolar Hamiltonians to the simpler XY 
(Planar) Hamiltonian^!?] ■ It is now clear, that even when one could control the 
many-body dynamics up to a high precision |13l IH| , its intrinsic chaotic insta- 
bility would easily degrade the quantum coherence^. This could explain the 
fast decay of the Magic EchoesfTB). the Polarization Echoes an d time 

reversal of cross-polarization. 19 appearing when one "inverts" the sign of the 
effective many-spin Hamiltonian. ^From this one learns that an attempt of dy- 
namical control of many spins systems should focus in the simplest interactions. 
Spins arranged in a linear chain with XY interaction have a clearly non chaotic 
dynamics. This would enable the transfer of the polarization amplitude from 
one edge to any specified nucleus in the chain. Besides, as the field of quantum 
computing with NMR |3U] grows steadily, one can foresee a new demand for 
alternatives to prepare and detect specific quantum states. 

The aim of this letter is to propose a novel concept that we call spin projec- 
tion chromatography. Given a collective spin excitation, it allows to filter out 
the components of a given total spin projection parity. To clarify the physics 
involved in this phenomenon, we use the known mapping between spins and 
fermions |^ E2] to develop a new formulation for spin dynamics based on the 
non-equilibrium formalism of Keldysh 23 . Its application to inhomogeneous 
chains and rings of spins 1/2 yields closed analytical expressions for the two 
site spin correlation functions that facilitate the prediction and understanding 
of interference phenomena in the spin dynamics at the molecular level. 



2 Spin Dynamics in the Keldysh Formalism 

The dynamics of a system with M spins under a Hamiltonian Ti. is usually 
described with the two site spin correlation function: 

Pf At) = ~ ■ (1) 

,,w (* eq \Si(t )Si(t )\* eq ) 

This gives the amount of the z component of the local polarization at time t 
on site /-th, provided that the system was, at time to, in its equilibrium state 
with a spin "up" added at i-th site. Here, Sj(t) = e lHt S^e~ lUt is the spin 

operator in the Heisenberg representation and l^eg) = 'YliN <lN i s the 

thermodynamical many-body equilibrium state constructed by adding states 
with different number N of spins up with appropriate statistical weights and 
random phases. The Jordan- Wigner (J-W) transformation |2] establishes the 
relation between spin 1/2 and fcrmion operators at site n, 



$n = ^ exp{i7r c+c TO }, 
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with c+ and c m , the creation and destruction operator for fermions, and 
the rising and lowering spin operator Si^ = S% ± iS%, where (u = x,y, z) 
represents the Cartesian spin operator. The relation 

c£ _ _ 1 

enables us to re- write the numerator of Eq. fl} as a density-density correlation: 

-l^(t)c f (t) - |c+(t )q(*o) + l* e9 > • (2) 

Since the last three terms in the rhs are constant occupation factors, the whole 
dynamics of the system is given by 



* e ,|c+(f)c/f)c+(M£;(to)|* e ,> (3) 



M („) 



E M 2 E(*t*° V-^(*)c/(*)^(*o)3i(*o) 

Af=0 7=1 w 



-0W 



Here, \E'^ V '^ and Z^r are the many-body eigenstates and the partition func- 
tion of the subspace with N particles respectively. We have assumed the high 
temperature limit, ksT = — » oo (i.e. e~^ w — > I), that leads to equal statis- 
tical weights |ajv | 2 = -pr (at) - Then, the invariance of the trace under cyclic 
permutations leads to: 

(*e,|a£(^/(*)#(*0)2i(t0) |*e,) = (*„e| ^ (*)?/(*) l*ne> • (4) 

The state |^ ne ) = G^C^o) \^eq) is a non- equilibrium many-body state generated 
at time to by creating an excitation at site i. The operator (t)cj (t) is evaluated 
over this non-equilibrium state. This is a particular case (t\ = t% = t and 
fc = I = /) of the particle density function^] in the Keldysh formalism: 

GZ tl fa,ti) = i (Vnelc+ihMh) \* ne ) . (5) 

The two-time particle density is related to the exact retarded (to. > ti) and 
advanced [ti < t\) Green's functions of the many-body system, 

G?,,(t 2 ,ti) = [G^(ti,* 2 )]* (6) 

= -i6[t2 - tl] (*ne|c fc (t2)^(ti) +^(ti)c fc (t2) |* ne ) , 

which for quadratic Hamiltonians are just density conserving single particle 
propagators. In more complex cases, they can decay through the many-body or 
environmental interactions. 



3 



Both Eqs. JSJ) and © can be split out into the contribution of each sub- 
space N. In particular, the non-equilibrium initial density at to describing an 
excitation at site i is, 

^fcj w) (*oj*o) = h {rrk^kj + jjE^ $k,i5i,i) , (7) 

a diagonal matrix involving equal weight for every site and an extra weight at 
the excited site i. In general, the initial state of Eq. (0) evolves under the 
Schrodinger equation expressed in the Danielewicz form |21 HI] , 

G <f ' (t, t) = 7i 2 £ ^ (*> *o)G < w (to, io)<T (to, t) 
l,k 

The first term in the rhs can be seen as an integral form of the (reduced) 
density matrix (p(t) = e -' m / h p{Q)e lUt / h ) projected over a basis of single par- 
ticle excitations in its real space representation as previously used in related 
problems |26l 1271 128| . However, the second term would collect incoherent rein- 
jections given by S < . These reinjection can compensate any eventual "leak" 
from the coherent evolution whenever the retarded and advanced propagators, 
G R and G A , contain decoherent processes. They also can account for interac- 
tions that do not conserve spin projection. Notice that Eq. JHJ) gives the freedom 
to operate with two time indices tk and ti, containing information about time 
correlations 00] - Hence, one could go beyond the standard initial value 

solution of the Schrodinger equation and include a continuous coherent injection 
of probability amplitude as is required to describe time dependent transport in 
mesoscopic and molecular devices jSOl EI]- Future applications may include the 
descriptions of the chemical kinetics at low rates of laser pumping into the ex- 
cited state[H21 Oni 12] and the gradual creation of homonuclear coherence in 
NMR cross-polarization transfer |5] . 

Equation (JSJ allows the evaluation of Eq. J3J which gives the non-trivial dy- 
namics of Eq. 0J. Using the high temperature limit, eq \c^ (t')c^(t') Y$ e q) — * 
|, the polarization is written as 

JV=1 

with 

pfi\t) = BG<f\t,t)-i. (10) 

Here, the dynamics PfP (t) of each subspace N, is weighted with the fraction 
of states that participate in the evolution. 
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We are now going to focus into systems where Eq. © takes the simplest 
form. If Ti commutes with the number operator, as is the case of the Heisen- 
berg and truncated dipolar Hamiltonians, one is sure that different subspaces 
N are decoupled. Further reduction could be obtained for Hamiltonians which 
are quadratic in the fermionic operator, such as the XY in proper conditions, 
as they might be reduced to non-interacting particles where S, , (iv) = 0. Such 
cases have fully coherent evolution and enable the best manifestations of quan- 
tum interferences. Besides, these one-body solutions are needed to build up a 
decoherent perturbative approximation of the many-body dynamics. 



3 Solution of Chains and Rings 

Consider a linear chain of M spins— i in an external magnetic field and inter- 
acting with their nearest neighbors through an XY coupling. Its Hamiltonian, 

M 

H cham = Ml n S^S n — ~ + \J n+ i. n {Sn + iS n +c.c). (11) 

n=l 

has a Zeeman part, 7i z , proportional to S* with fl n the chemically shifted 
precession frequency; and a flip-flop term, H-xy^ where J nn+1 is the coupling 
between sites n and n + 1 with J M 1 = 0. 

Due to the short range interaction, after a J-W transformation the only non- 

C n+l C n 

\ C'4-a4-!-iC' fif omn nmio/iliVin / \ 1 Q ^ 

■ N 

non-interacting fermions. The eigenfunctions 



zero coupling terms are proportional to c^ +1 c„ = S^ +1 S n . Each subspace with 
(^f) states of spin projection (j2n=i ^n) = N — M/2 is now a subspace with N 

\E , y V '*\ are expressed as a single 



Slater determinant built-up with the single particle wave functions (p a of energy 
e a - Under this condition, G^ N) (t) = Gf^t) for every N, and Eq. reduces 
to 

Pf*T(t) = h 2 \Gf ti (i)| 2 . (12) 

Even for inhomogeneous systems, this can be solved analytically or numerically 
at low computational cost. 

For spins arranged in a ring, the closure of the chain requires an extra 
term in the Hamiltonian of Eq. ^] preventing the use of Eq. After the 

J-W transformation this term takes the form , 

M 

\Jm,\S^S m = -| Ji, M ctc M exp[i?r c£<J m ] . (13) 

m— 1 

Thus the simple "flip-flop" interaction of the spin system becomes a "complex 
many-body" interaction term in the fermionic representation. The exponential 



can take the values ±1 depending on the parity of N = (\& 
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N ^ N) ), the 
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TIME [2fi/ J] 

Figure 1 : a) Evolution of the polarization at each site of a ring of size M = 20 af- 
ter a local excitation at site 10. b) Contour plot of the polarization. Black filling 
represents P > 0.09; white is P < 0.045 and grey corresponds to intermediate 
values. 
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integer number of particles. Thus, there are two possible one-body Hamilto- 
nians determining the wave functions ip a (n) used to build the Slater determi- 
nant tyy. All subspaces with the same occupation parity are described with 
the same Hamiltonian. The different signs can be seen as a parity dependent 
boundary conditions (i.e. either periodic or antiperiodic) or even better as the 
effect of an effective potential vector A, expressed in a singular gauge: 



M 

V c+c 



od2 



Adl 



4> 



i- <t Adl =i2?r 

hj $ 



(14) 



Here, the second equality defines the singular gauge for A which modifies the 
hopping term according to the Peierls substitution j85j . This results in an 
Aharonov-Bohm phase 27rJ^, where $ is the "fictitious" magnetic flux pierc- 
ing the ring with value $ = or $ = i$ = \hje depending on the particle 
number (spin projection) parity. Therefore, an odd occupation (N = 2n + 1) 
yields $ = 0, while the even occupation (N — 2n with n integer) corresponds 
to $ = $0/2. The resulting Hamiltonian, expressed in the singular gauge and 
satisfying c n = c n+M , reads: 

H ring ($) = H chain + ±( Jim^'^Cm + cc). 
The resulting propagators, 



Gf/ n+1 \t) = G% (i;* = 0) and 



G 



R (2n) f ^ _ 



(t) = G* (t; i$ ), 



are inserted in Eq. J5J) to obtain the polarization from Eq. J5J: 

jring/^ _ ft 2 \\ r ,R /. ^ ~nI 2 



(15) 

(16) 
(17) 

(18) 



Again, the dynamics is completely solved by evaluating the single particle 
Green's functions. 



4 Spin Projection and Excitation Velocities 

Let us focus into the spin dynamics of homogeneous rings (J n +x,n = J and 
Q n = 0). All the properties can be expressed in terms of the single particle 
wave-functions which we solve in the Coulomb gauge A — x r obtaining 
<Pk{n) — (0| c n c^ |0) = exp[ifen]/\/M; here |0) is the particle vacuum and 
creates a plane wave with a kinetic momentum hk — hk + <f^^gr obtained 
from the stantard discrete wave vector given by k — s2tt/M with s = 1, . . . M; 
the lattice constant is a = 1. We recall that $ is a fictitious flux identified 
with the parity of the number of spins up in the original problem. Hence, a 
high temperature excitation will require solutions with both fluxes $ = and 
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Figure (JIJ shows the evolution of a local excitation through every site of a 
ring with M = 20 sites at high temperature. The original excitation at the 10-th 
site G< m (t Q ,t Q ) — j^8 n ,w&w,m splits-out into two wave packets with opposite 

velocities (i.e. + or -). Formally, G< m {t,t) = Gnjm(t,t) + G^m [t, t) where 
the partial particle densities propagate according to the definite sign in the 
kinetic momentum hk. Let's call j n +i,n(i) = j„ \t)— j^+iit) the current density 
between sites n and n+1. From Eq. (5.1) in Ref. [251 we § e ^ the partial currents: 

Jn^Ht) = 5 ^«,n+l^n±l, «(*>*) - \j n+ i^ n G^^ ±1 {t,t), (19) 

defining the propagation in each direction. The components with maximum 
velocity become increasingly dominant in the determination of the front of the 
wave packet. The eigenenergies er = Jcos[fc] allow the evaluation of the group 
velocity v~ — jr J| . The partial current and the corresponding local density can 
be used to evaluate the mean velocity as v = j/p, i.e.: 



Ln T Lr ">' 1 \ l i X, J 



0<k 2 <ir/a 

setting a lower bound for the peak velocity. The revival time Tme for the 
mesoscopic echo is proportional to the ring size M and inversely proportional to 
v and can be estimated as M/i> max < Tme < M/v. In Fig. 0Jj) the dominant 
group velocity manifests as a linear dependence of the wave packet maxima with 
time. 

There is a fundamental difference between rings of odd and even number 
of sites M. For odd sizes, the eigenstates determine the same set of velocities 
regardless of the boundary condition determined by the Aharonov-Bohm flux 
(which is half multiple of $ = $o/2 or $ = 0). i.e. neither the available energies 
nor the velocities depend on the parity of the particle number, see the top of 
Fig. J2J). However, the bottom panel shows that for even sizes, this symmetry 
is always absent and the set of velocities depend on i.e. on the parity of the 
particle number. 

Different velocities for the eigenstates lead to different mean velocities for 
the localized excitation. Indeed, the evaluation of Eq. (|20f) for even values of 
M yields 

, $ 2tt n J 2 , $ 2?r x 
v($) = w(0)cos( — — ) - TTrsin — — , 

while for odd M gives a mean velocity which does not depend on $ (i.e. on the 
particle parity) 

v($)=v(0), with ^ ) = -^A/tan(^). 

It is evident that the mean velocity for the even ring sizes satisfies v(Q) > 
v(0) indicating a faster propagation of the excitations with an even number 
of particles (i.e. in the subspace with even spin projection). For M = 8, 
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Figure 2: Dependence of velocities of eigenstates for rings with odd (7) and 
even (8) number of sites as function of the magnetic flux $. Full and empty 
dots are eigenstates for the cases $ = and $ = $o/2 which corresponds to 
odd and even total spin projection parity respectively. In rings with odd spin 
number parity the velocities do not depend on spin projection parity. 
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t>(i$o)/f(0) = 1.08, in contrast, for M = 7 this ratio is 1. This result is 
illustrated in Figure © for an M = 8 spins system . Panel a) shows the even 
^(evcn) _ ^ ^(2n) | gz |^(2n)^ and odd components of the local polarization 
for the ring with eight sites. Panel b) shows the total polarization represented 
by the sum of both contributions. In summary, in rings with even number of 
spins an excitation in a subspace with even spin projection moves faster than 
one in a subspace with odd spin projection. This conclusion, easily derived in 
the language of "particles" and "fluxes" , would have been more elusive in the 
original language of spins. The novel concept revealed by the previous simple 
example is the possibility to manipulate collective entities as those constituted 
by the allowed total spin projections in a many-spin system. 




Figure 3: Spin Projection Chromatography in a ring with M = 8. Panel a) 
shows the contribution of even and odd spin projections to the local polarization, 
the sum of both, is shown in panel b). 
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5 Conclusions 



We must emphasize that the transparency of the formalism fully presented here 
for the first time, Eqs. 15l!)l) . has already been instrumental to reveal the simple 
concepts hidden in the various problems related to spin dynamics such as the 
nature of Mesoscopic Echoes |51 I10| and Loschmidt Echoes |15) . 

The experimental realization of an XY Hamiltonian can be achieved through 
the truncation of the more complicated Heisenberg (J-coupling) Hamiltonian 
by the pulse sequence developed in Ref. \T2\. Its application to relatively 
small cyclic molecules of even parity (e.g. a derivative of octatetraene) would 
enable the observation of the described dynamics. The evident difference in 
the delay time of the ME of each spin projection allows one to manipulate 
them independently. One could conceive filtering out components of the total 
spin projection of a given parity. Consider a local magnetization (aligned in 
the z direction) at a site distinguished by its chemical shift in the M = 8 ring. 
This can be achieved by selective irradiation of one spin or by other polarization 
transfer method. This non-equilibrium initial state starts an evolution described 
by Eq. (|18|) . After a time Tme, the excitation will return to the initial site but 
the components moving in an even spin projection space will come back first. It 
is then possible to give a selective 7r/2-pulse (purification pulse) at t = 7.5(2h/ J), 
when the even spin projection has its maximum probability amplitude over the 
site and the corresponding probability amplitude for the odd spin projection 
is still close to cero, see Fig. (|2Jl. Thus, the even total spin projection has 
been tumbled down towards the xy— plane by the purification pulse leaving 
the odd one unaltered, i.e. an almost pure odd total spin projection state lies 
along the z— direction. Hence, the collective nature of the excitation can be 
used to achieve a dynamical filtering of the parity of the total spin allowing the 
purification of the desired component. This justifies the name of spin projection 
chromatography. 

We acknowledge the financial support from ANPCyT, SeCyT-UNC, CON- 
ICET and Fundacion Antorchas. 
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